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The dispersion relation is derived for the coherent waves in fiuid or elastic media supporting 
viscous and thermal effects and containing randomly distributed spherical scatterers. The formula 
obtained is the generalization of Lloyd and Berry's [Proc. Phys. Soc. Lond. 91, 678-688, 1067], the 
latter being limited to fluid host media, and it is the three-dimensional counterpart of that derived 
by Conoir and Norris [Wave Motion 47, 183-197, 2010] for cylindrical scatterers in an elastic host 
medium. 

I. INTRODUCTION 

This study considers wave propagation through a homogeneous isotropic host medium containing a large number 
of randomly and uniformly located spherical scatterers. It is well known that for small enough concentrations the 
physical medium may be replaced by an effective homogeneous medium in which coherent waves propagate. These 
waves correspond to the average, over all possible locations of the scatterers, of the multiply-scattered field in the 
actual host medium. Our concern here is to find the effective wavenumbers of the coherent waves at low frequency 
and low concentration of spheres. This was achieved some time ago by Lloyd and Berryi for the case of acoustic 
waves with P = 1 (P: number of waves that propagate in the host medium). A new and clear derivation of Lloyd 
and Berry's formula was recently given by Linton and Martin^ using a procedure developed earlier for cylindrical 
scatterers in Ref4. They considered a uniform concentration of scatterers satisfying the hole correction of Fikioris 
and Waterman^, and found solutions of the implicit dispersion relation as an expansion in terms of the concentration 
of scatterers, up to order 2, under the low frequency assumption. Generalization of Linton and Martin's formula for 
cylinders in a fluid^ has been given by Conoir and Norris^ for cylinders in an elastic solid with P = 2 (compressional 
and shear waves propagate in elastic solids). The present paper is an extension of the results of Linton and Martin^ 
to homogeneous isotropic host media supporting viscous and thermal (damped) waves for which P = 3. Examples of 
such media include viscous fiuids^ii with compressional and transverse waves (P = 2), and thermoelastic solids2r2S 
in which compressional, transverse and thermal waves propagate (P = 3). The results obtained here can serve as 
a starting point to generalize effective medium theories that include multiple wave interactions, a broad frequency 
range, and finite concentration levels, such as the ECAH model (Epstein, Carhart, Allegra and Hawley)^. 

The analysis begins in section |TT] with a description of the scalar potentials that are used to describe the wave 
propagation in the host medium. The choice of the potentials depends on conditions of symmetries which are discussed. 
Multiple scattering theory is used in section IIIII to derive the modal equations (Lorentz-Lorenz laws) relating the 
amplitudes of the P-wave types. In section IIVI the modal equations are written in matrix form which allows us to 
obtain a compact form of the dispersion equation governing the coherent waves. The low concentration and low 
frequency assumptions are then introduced in section |V] and expansions of the effective wavenumbers up to second 
order in concentration are derived. The effective wavenumbers are expressed in terms of series, and, for the faster 
wave, acoustic wave in viscous fiuids and compressional wave in viscoelastic solids, the wave number is given by an 
integral relation that generalizes the one of Lloyd and Berryi. 



II. THE DIFFERENT TYPES OF COHERENT WAVES 



We first define the multiple scattering problem for identical spheres in homogeneous isotropic host media in which 
three type of waves are present, compressional (c), viscous or shear (s) and thermal (th). The multiple scattering 
theory used here follows the lines of that first developed by Fikioris & Waterman^ for acoustic waves. 
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A. The different types of waves in the host medium 

The dynamic displacement field u in the host medium may always be decomposed as w = VtAl + V x -0^, with 
the scalar potential for longitudinal (compressional or thermal) wave motion, and ij^R the vector potential associated 
with rotational waves^. Each wave type propagates with its own complex wavenumber kp, with p an integer numbering 
the wave (1 < p < P ). The rotational wave field itself can generally be partitioned into a shear and a transverse 
wave, corresponding to the decomposition of the vector potential ipR as the sum of two orthogonal vectors V'i?, = 
V X (ipre^r) + (1/fcp) V X V X (xre'r) with the unit vector along the radial coordinates^^. The Debye potentials ip 
and X s-i'S associated with the shear (s) and transverse (T) waves, respectively. As a consequence, as many as four 
distinct wave types can propagate in the host medium (p = c, th, s, T). 

Our objective here is to determine the dispersion relations of the coherent waves in the presence of spheres that are 
randomly and uniformly located in the host medium. A priori^ there should be one coherent wave for each longitudinal 
wave, and one for each shear and each transverse wave. This is a natural hypothesis if the scatterers are scarce, that 
is of low volumetric concentration. It is important to note that the dispersion equations are a characteristic of the 
effective medium itself, and, as such, they should not depend on the type of the incident plane wave, which will 
be supposed from now on to be a compressional wave. This choice implies symmetries on the scattering by one 
scattereri^iiSj which, in turn, have implications on the form of the analytical expressions for the potentials associated 
with the coherent waves propagating in the effective medium. 

B. Symmetries 

The scatterers are assumed to be randomly and uniformly located in the semi-infinite region z > oi the host 
medium, with an harmonic compressional plane wave incident on the boundary z = 0. Coherent wavenumbers do not 
depend on the angle of incidence of the incident wave^, and so with no loss in generality we assume normal incidence. 
Let us first consider scattering by a single sphere of an incident compressional plane wave propagating in the z direction. 
Due to the symmetry of the sphere, the scattered fields do not depend on the azimuthal angle ip, with spherical 
coordinates defined by the sphere center {(x,y,z) = r {sin cos (p, sin 9 sin (p, cos 0))^^. Consequently, the scattered 
transverse T wave with pure azimuthal displacement is zero. The question now is whether a transverse coherent wave 
can arise from the multiple scattering process. The answer is not straightforward, because the independence of a field 
with respect to the azimuthal angle in a given spherical coordinate system does not guarantee independence of that 
same field with respect to the azimuthal angle in another spherical coordinate system. This means that each multiply 
scattered field is a combination of longitudinal, shear and transverse waves. However, if the scatterers are uniformly 
distributed, the average introduces an homogenization of the different fields in the medium, and the effective fields 
exciting a given sphere should not depend on the azimuthal angle ip, and, thus, the coherent transverse T field should 
vanish. 

Another way of considering the symmetry is to bear in mind that coherent waves result mainly from constructive 
interferences of waves traveling from one sphere to another along a "straight line" (the direction of propagation of 
the coherent waves). The effective coherent fields maintain the same symmetry as that for the field scattered by a 
single sphere. We therefore conclude that the coherent fields do not depend on the azimuthal direction cp in the {x, y) 
plane perpendicular to the direction of the incident plane wave. Mathematically, this result is exact up to order 2 in 
concentration, as shown in Appendix A. In the following, the longitudinal fields are described by scalar potentials of 
displacement, and the shear fields by the only non-zero component (azimuthal component) of the vectorial potential 
of displacement related to shear waves. These scalar quantities are denoted by (z?^^-* , with p denoting again the type 
of the wave {p = c,s, th). 

III. MULTIPLE SCATTERING EQUATIONS 

Harmonic wave motion is considered with time dependence exp(— iwt) understood. The notation c = 1, s = 2 
and th — 3 is employed, similar to that used in Rcf^. The potential function ip^g\r; fj) represents the wave of type 
p scattered by a target centered at r) and observed at r; the potential ip'^^\r;fj) denotes the field of type p that 
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excites a scatterer centered at Tj and observed at r. The following fundamental identity, which is a straightforward 
generalization of Eqs. (6) from Refi^, provides the integral equation governing the coherent fields, denoted by brackets, 
in the presence of a uniform and random array of identical scatterers: 

(^(r)(r;rl))= <5,i^W(^1+E / Ar,n{r-,r^)T'^^{r,){^^i\r-r,)). (1) 

The integration in Eq. ([T]) is over the semi-infinite region {z > 0) containing the spherical scatterers. The function 
n{fj, ri) is the conditional number density of spheres at rj if one is known to be at ri, see Ref.— . In the following we 
assume a constant density no of scatterers of radius a, and conditional number density given by the hole correction 

\ nn for |r — r, |>6, with b > 2a. 
^ |0 otherwise. ^ ' 

The incident comprcssional wave is assumed to be a damped plane wave 

^»c(^i = e'^''^ with kp = fc; + tk'; (fc; > o). (s) 

The potential functions are expressed as infinite series of spherical harmonics that respect the underlying symmetry. 
For the reasons described before they do not depend on the azimuthal direction ip, and so 



{^E\r,rj)) = E Alr\rj)jnikpPj)Pn{coaeip,)) with p, ee r - r,. (4) 

ri=0 



It proves useful to introduce shorthand notation for spherical harmonics (not to be confused with ordinary Bessel and 
Hankel function notation) 

Jn{kp,Pj) ^ jn{kpPj)Pn{cOs9{pj)), Hn{kp, Pj) = h^^^kpP j) Pn{ COS 9 (pj)) , (5) 

plus the following definition for the action of the transition operators T'''^ on a spherical harmonio^iii^ 

T^"" (rS- ) J„ (fcp , ) - if„ (fcp ,/?,). (6) 

The modal coefficients T^p in Eq. ^ can be numerically calculated following the procedures of Refsi^i^. Equation 
dD) can now be expressed, using Eqs. ^ and (jH), 

+ 00 +CJO P 

Y,Alr\ri)Jn{kp,pi) = Ve*'^^^ +EE / dr,n{r-,rMi'\r,)TrHn{kp,p,). (7) 

Decomposition of the damped incident plane wave into spherical harmonics in the coordinate system centered at fi 
leads to 

+00 

[Ai^Hn] - 6pii"{2n + l)e'''^^-^] J„(fcp, p,) 

n=0 

+00 P 



Y^Y. I <irjn{r-,fMl'\r,)T:iPH„ikp,p,). (8) 

n 1 ^ 



ri=0 q=l 



The addition theoremi^ allows us to write the series in the right hand side of Eq. (|5]) as a function of the coordinates 
centered on fi. 



+00 +z> +00 

ffn(fcp,p,) = E E E(-i)''*''"'""(2^ + i)G(0'^;-/^''^;^)^ 

xhl^\kprij)p-^'{cose{ri,))j,{kpPi)Pli{cos9{pi)) with ry=rl-rj, (9) 
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where the Gaunt coefficients G(0, n; —fi, v] m) are defined by 

oo 

P^' { COS 0) { COS 0) = ^G(m,n;Ai, I/; €)P;"+''( COS 61). (10) 

e=o 

By assumption, the incident damped plane wave impinges on the (z = 0) interface at normal incidence, and gives 
rise to damped coherent waves that propagate and are attenuated in the same direction z. Accordingly, we seek 
coherent plane wave solutions obeying the Snell-Descartes laws of refraction. We thus search for the solutions of Eq. 
([5]) in the form 

p 

M?Hrj) = z"(2n+ l)^^(«P)e'«''^^. (11) 
p=i 

Substituting from Eqs. ([9]) and pTjl into Eq. (|8]), making the change of variables of integration fj = fx — rij, and 
noting that the integration over clipij ~ d(p{r\j) with drij ~ sm0ijdrijd9ijdipij gives rise to zero except for fi ~ 
because of the term exp[ — i^(^(rij)] , yields 

P -fCJO +CXD 

= j^^e^'^''^^ +^^^i^(2i.+ l)TJP4«'=)G(0,t.;0,n;£)4^^ (12) 
5=1 j/=o e=o 

where the quantity I^'''^ is 

ll^HCk)^ J dr,n{r-,n)e'^''^h['\kpn,)Pe{cos0irrj)). (13) 
Taking into account the hole correction ([2]), and following the analysis of Ref.— , we obtain 

with (14a) 



(p) _ 2noTTi 



26 ,jn) , . ^ ic. i 

re" 



N^/\a) = abji{^kb)h[^\kpb) - kpbMab)h['^' (kpb). (14b) 

Inserting Eq. ([T4| into Eq. (fT2|). and equating the coefficients of exp(i^fcZi) to zero gives what is known as the 
Lorentz-Lorenz law (equating the coefficients of exp(ifcpZi) to zero gives the extinction theorem). Finally, we obtain 
the equations for the amplitudes yli^*^'' , 

. , P +00 +C30 

^(,,) ^ 1^ ^ ^ E(-l)'(2- + l)Tr4'^)<(a)G(0, 0, n; I) (15) 

^fe P q=l u=0 (=0 

with 1 p P and 1 A: P. 



IV. MATRIX FORM OF THE MODAL EQUATION 

We now seek to write Eq. pSI) in matrix form in the same way as in Refj^ii^. With no loss of generality, consider a 
given effective wave number {k — c, s,th). Equation (jlSp can therefore be written in simplified form by dropping 
the index k (S,k = S,)- Define the non-dimensional parameters 

?/p = ^■^ — fcp and e = — 4mo , (16) 

and the use of the relationi^ 

(-l)^G(O, i^; 0, n; i) ^ (-1)'^+"G(0, i^; 0, n; i) (17) 

implies (1 ^ p ^ P ) 

, P +00 +00 

^(p) _ !£!^^^^(_l)-+»(2^+ l)T74«)iV(^)(OG(0,i/;0,n;€) = 0. (18) 
yp o=i u=o e=o 
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In order to reach our objective, we define the infinitely long vectors |e), (e| = |e)* and square matrices T'^, q''^^ from 
their components, e„ ~ (^1)" a-nd 



+00 



TZ = 5n.(2i^ + 1)TJ^ g(fj(e) = J^[^k,bY,N^P\OG{0,^y;0,n;i) - l](|e)(e|). 

dp) 



(19) 



We also introduce the unknown vectors jop) with components An , the vectors \ep) consisting in the combination of P 
zero vectors except at the p*^ place where their components are equal to yjT^kp\e), the |a) vector that is a collection 
of the jop) vectors, and the block matrices 





I 





0" 




Ip = 





I 





T = 










I 





rj-ill r|-i21 r-|-i31 
r|-il2 r|-i22 r-|-i32 
rpilS rpi23 fTi33 



Q 



Q 



(1) 
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(3) 



q'^^ 



q: 



(20) 



Using these definitions the system of equations ([T5|) can be replaced by the equivalent condition 



{ip-eQT-e^ 



|ep)(epl 
Vp 



T}|a) = |0). 



(21) 



The matrices and vectors as defined are based on the assumption P — i. For smaller values of P, or equivalently, no 
coupling between the wave types (T'^^ = T''^(5qp), Eq. (I^Tt decouples into separate equations for each wave type, each 
analogous to the acoustic (P = 1) case. 

Equation (j2ip embodies the multiple scattering of the P wave types in a single consistency relation. Its structure 
has been specifically chosen so that it is the same for spherical and cylindrical geometries^ whatever the number of 
waves (P ~ 1,2,3). We can therefore follow the procedures as developed in Refi^ii^ for the cylindrical case. After 
multiplication on the left hand side by and introduction of the following quantities 



we get 



|/p)=T^|ep), (/p| = (ep|T^ \h) ^T--\a) and Q = T^QT^ 



{Ip-e(Ip-6Q)-if]MM}|5)^|0). 



Vp 



(22) 



(23) 



The modal equation is obtained by setting to zero the determinant associated with the infinite homogeneous system 
of equations Eq. (|21[> ^i^. This can be written in a much simpler form by writing Eq. (|23p in the form 



{Ip + ^|5p)(/p|}|6) = |0), where |gp) ^ -£y^-i(Ip _ eQ)-i|/p). 

p=i 



(24) 



Noting that the matrix in (j24p is the sum of the infinite dimensional identity plus a matrix of rank P, the determinant 
can therefore be reduced to one for a P-dimensional matrix through the use of the identity 



det I Ip + ^ l5p)(/p| ) ~ det (ip + m) where 



ma 



''qp ~ {9q\fp)j 

and ip is the identity matrix of dimension P. Setting the determinant in Eq. (j25p to zero yields 



yi 


1/2 


2/3 


_gMl2 

VI 


V2 


^ A/32 

2/3 


VI 


^ A/23 
V2 


ys 



= 0, 



(25) 



(26) 
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or equivalently 



where the matrix elements are 



yi - eMii -eA/21 -eMai 
-eMi2 2/2 - e^^22 -£^32 
-eMi3 -eM23 2/3 - e^^33 



M,p(0-(/J(Ip-eQ)-Vp) 

+ 00 



(27) 



(e,|T|ep) + 5](e,|(TQT)"|ep)e'^ 

+00 



(28) 



n=l 



This indicates that the elements of the modal equation Eq. (P7)) can be calculated without evaluation of the square 
root matrix T^. It is also evident that the above reduction in the size of the system determinant is valid as long as 
the matrix Ip — eQ is non-singular, which is certainly valid for small e. Equation (|27p is the fundamental equation 
for determining the coherent wavenumbers ^p. 

V. ASYMPTOTIC SOLUTIONS OF THE WAVENUMBER EQUATION 

This section considers asymptotic expansions of the solutions, valid in different limits; first for low concentration 
and then for low frequency. 



A. Low concentration expansion 

Rather than working with the wavenumber directly it is more convenient to expand the solutions of Eq. (j27p 
about one of the three leading order solutions j/p = ( ^ = fcj,). The non-dimensional parameters e is small at low 
concentration, |e| 1, and we therefore assume a formal asymptotic expansion in e: 



Inserting the asymptotic expansion into the modal equation (|27p provides (c/. Appendix B) 

yW=M(°)(fcp), 



;f =M«(M + E 



Mi°^(fcp)Mj°^(fcp) 

1,2 _ 1,2 

ftp 



It follows from Eqs. ([TH) and (HH) that 



m(0) = ^ 



E(2n + l)Tr, 



3 00 00 



(29) 

(30a) 
(30b) 



(31a) 
(31b) 



q=l n=0 u=0 
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with 

+00 



X [fcp6j;(fcp6)/if - kgbMkpb)h^l'>' (kgb)] L g ^ p, (32a) 



+ CXD 



e=o 



+ Yb " + 1)] J^(fcp^^)/if ^(fcp&)| (32b) 



where the following relation has been taken into account, 

+ CXD 



Y,GiO,i^;0,n;e)^l. (33) 



B. Low frequency expansion 

The long wavelength limit is defined as kpb — > 0, in which case Eqs. ((32|) reduce to (see Eq. (p3)) ) 

+00 , 

Q^r?} ikp) = , " pJ -ir'^ E GiO, 0, r.; ^) 1] , g ^ P, (34a) 

+ 00 

^-(M = :^(-l)"+''E^G(0''^;0'";^)- (34b) 

The low frequency version of Eq. (I31b[) is 

2 +00 +00 +00 

M^lHkp) = ^ E E E(2" + + 1) ^(0, 0, n; £) 

P n=0 i/=0 e=o 



X 



l^rr^r + E k (I? k^) fc")' - (35) 



9#P 

and hence the low frequency expansion of the effective wavenumbers are 

>2 00 o 2 2 +°° +°° +°° 

=1 - E + - ^ E E E(2" + 1)(2- + 1) G(0, 0, n; ^) 

P P n=0 P ri=0 !v=0 i=0 



iT^PT^p + Y^ 2fcp (^)^rrr,^. ^. (36) 



If the host medium is an ideal fluid, i.e. P = 1 as in Ref3., Eq. simplifies to 

f2 oo 

|=l-4-SE(2- + l)^^r 



1,2 1,3 

o 2 2 +°o +°o +°o 

- ^ E E E(2" + 1) + 1)^(0' 0' ^) (37) 

P ji=0 i/=0 £=0 

which is exactly the same relation as that obtained by combining Eqs. (1.2,4.29, 4.33,C.4) from Refi^. Equation 
pop contains additional terms which are not in Eq. ([37]) and are clearly connected to the coupling between the 
compressional, shear and thermal waves {q =/= p). These coupling terms are neglected in the EC AH model (Epstein, 
Carhart, Allegra and Hawley)^^. Equation can therefore serve as the starting point for further developments, 
such as the calculation of the Rayleigh limit, for example, and is the principal result of the paper. 
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C. Generalization of the Lloyd and Berry formula 



The advantage of the Lloyd and Berry formula is that it expresses the wavenumber at low concentration in terms 
of the far-field scattering function only, rather than the T-matrix elements. In the present context, this requires that 
we express series in (j36|) as integrals of the far-field scattering functions defined as follows 



+00 



= Y.i2n + l)T^Pp„{cose). (38) 



n=0 



The analogous acoustic problem involves only the terms in Eq. p7p which have been shown to be equivalent to the 
following expansion in the concentration^ 

= + Smo + S2nl, (39) 

with 

4-7 TT 

Si = -^/^^(O), (40a) 

^2 = ^1 [F^of - [F^{7:f+ £ (40b) 

where S2 in Eq. (j40b[) is the formula initially given by Lloyd and Berrjii. In order to express the general coupled wave 
problem in a similar form it is necessary to represent the final series in Eq. (j36p {q ^ p) as integrals of the far-field 
scattering functions. 

The relevant series that appears in the multi-wave system is (see Eq. (I36|)) 



+00 +00 +00 J 
^(k) = ^^^(2n + l)(2i. + l)T7rP««;^G(0,i/;0,n;£) where ^ = (41) 



n=0 i/=0 £=0 



Two cases need to be considered depending upon whether |k| < 1 or |k| > 1. For the case in which |k| < 1, we can 
define the function 

-|-oo 

g{K,9) = K'"(2m + l)P,„(cos6'). (42) 



The use of the following two relations (c/. RefJ^) 



Tn=0 

15^ 



2 

F„(cos6')P™(cos6') sin6'd0 = -(5„m, (43a) 

2?! -I- 1 



-t-00 



Y,^''Pn{coa0)= ^ -r with \k\<1, (43b) 

(1 — 2KC0St' + K^)2 



then provides, respectively. 



.g(K,6')P„(cos6l)sine'd6l, (44a) 



g(^,g)^ [1 + 2^-f] \ -r. (44b) 

dK' [I -2k cos e + K^) 2 

The product of the far-field functions follows from the definition (|38p as 

+00 +00 +00 

/«P(0)/^"?(0) = ^^^(2n-M)(2:. + l)T,frrP£(cos0)G(O,i^;O,n;^). (45) 

n=0 u=0 1=0 

Combing the above results and performing the differentiation in ()44bp gives 

S{k)^1-{1-k^) f f''P{e)fP''{9){l-2KCOse + K'^yhme de. (46) 

2 .In 
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The series S{k) is therefore convergent since it can be expressed as the integration of a continuous function over the 
compact interval [0, tt]. Although we did not manage to express the remaining series in Eq. (j36p in terms of the form 
functions when |k| > 1, the series are always convergent because there is only a limited number of significant values 
of the T^P coefficients, see the discussion on this point in Rcf.^. It is also easy to show that the Gaunt coefficients 
decrease much more quickly that the function increases with £ +oo even for |k| > 1. 

If p indicates the faster wave propagating in the medium, then it is always true that |k| = \kp/kq\ < 1 whatever the 
value oi q p. This is the case for acoustic waves in viscous fluids and compressional waves in viscoelastic solids. So, 
in such cases, if is the wavenumber of the acoustic or compressional wave, it follows in a straightforward manner 
from Eqs. (jMl) to gG]) that 

^p = kl + Sina + S2nl + S^^'^^nl (47) 

with 

q^p'^ {kl + kl-2kpkq COS 6)^ 

The additional term 5^2^ indicates the coupling between the compressional (p), shear and thermal waves {q ^ p)- The 
formula given by Eqs. (|47p combined with (j48p generalizes the identity for acoustic waves derived by Lloyd and Berry^ 
for which (Sj — 0. It is of interest to compare the structure of ^2 with that of the corresponding term for the effective 
quasi- longitudinal wave in the presence of cylinders in an elastic host medium^ (Theorem 1). The coefficient for the 
latter case is obtained from (j48]) by removing smO in the numerator and replacing the power | in the denominator 
with 1, that is, by making the obvious changes one would expect for 2D as compared with 3D. 



Appendix A: Azimuthal waves at second order 

The aim of this appendix is to show that the coherent fields do not depend on the azimuthal direction tp up to 
second order in concentration. We refer here to equations in the paper of Linton and Martin^. The fact that it only 
considers compressional or acoustic waves [P = 1) is not important for the present demonstration, which is the same 
for one or for several waves {P = 3). We begin with the modal equation for oblique incidence, (c/j^ (Eq. 4.20)) 



»no(4^)2(-l) 
/fc(fc2 - K^) 

with the spherical harmonics and the Gaunt coefficients defined by, respectively. 



^"+ '"fc j- j^2 E E E^-^'^''^r™(^)N,(7^6)G(n,m;z.,-M;9)=0, (Al) 



= YrTie.^) = ^-^r\j'^^\I^^^^^Pn\^ose)e^"^^, (A2a) 

oo 

yr(r)y,^(?) ^ n"^^(?)G(n, m; fi; q). (A2b) 

q=0 

Gonsequently, it is necessary to prove that the unknown coefficients are zero except for m = 0, at least up to 
order 2 in concentration. We use the following result for the expansion of up to and including the second order 
in concentration {cf^ (Eq. 4.20) and its preceding equation) 



F™=F:(K)F + no{y::(K)V^+i^(-l)™i^^ J2 Z,Y':{K)XZ'} (A3) 

1/—0 fj,— — i/ 

with^ (Eq. 4.27) 

+ 00 

XTJ" = E 5T""(K)G(n, m; -p.- q) dq{kb), (A4) 

q=0 
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where F is a constant and where V— (Eq. 4.31) and dq{kb)^ (Eq. 4.23) are functions which do not depend on the 
indices n and m. Since the coherent wavenumbers do not depend upon the angle of incidence*^, we may consider the 
case of normal incidence for which K = (0, 0, 1)^ (Eq. 3.15), so that 



2q+l 

( 

An 



Equation (jA3p therefore reduces to 



+ 00 +Z.' +00 



^ Z.^/i2v + l)(2g + 1)5^0 S^m G{n, 0; i^, 0; q)dq{kb), 

v—O fj,— — !^ q—0 



which proves that = if m 7^ up to second order in concentration. 

Consequently, at normal incidence and up to order 2 in concentration, Eq. (|A1[) becomes 



F' 



»"-o(47r) 



2 +00 -\-oo 



.o./2g+l 
47r 



N,(/a)G(n,0;i/,0;g) = 0. 



(A5) 



(A6) 



(A7) 



Our results correspond exactly to those of Ref»^. Indeed, we can easily show that Eq. (jA7p and Eq. ([T8|) with P = 1 
are exactly the same if we make the following identifications: k = ki, K = = Zi, = —T^^, ISIqlKb) = ikibNn\s^) 
and 



(A8) 



Appendix B: Some asymptotic expansions 



With no loss of generality we consider the first root yi. It follows from Eqs. (|16l) . (P^l) . and (PH)) that the asymptotic 
expansion of the 1 1 element of the matrix in Eq. (P?]) is 



y, - eMn = [y['^ ~ Mi°\kp)]e + [yf ^ - Mi\\kp)y + ... . 
Hence, to leading order Eq. ([27]) becomes 



(Bl) 



yW-Afi7(fcp) 
—M^2(kp) ki — 



-M[°\kp) kj~kl 

which implies the identity pOap . Inserting the latter into Eq. (jBip and Eq. ([27| then gives, at the leading order, 

y^^^~M[l\kp) ~A4fikp) ~A4fikp) 



= 0, 



(B2) 



^12 i^p) ^'1 ^2 





''■1 ''■3 



(B3) 



from which Eq. (j30bp follows. 
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